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ABSTRACT 

Using galaxy group/cluster catalogs created from the Sloan Digital Sky Survey Data 
Release 7, we examine in detail the specific star formation rate (SSFR) distribution 
of satellite galaxies and its dependence on stellar mass, host halo mass, and halo- 
centric radius. All galaxies, regardless of central-satellite designation, exhibit a similar 



bimodal SSFR distribution, with a strong break at SSFR « 10 



yr 



and the same 



high SSFR peak; in no regime is there ever an excess of galaxies in the 'green valley'. 
Satellite galaxies are simply more likely to lie on the quenched ('red sequence') side of 
the SSFR distribution. Furthermore, the satellite quenched fraction excess above the 
field galaxy value is nearly independent of galaxy stellar mass. An enhanced quenched 
fraction for satellites persists in groups with halo masses down to 3 x 10^^ Mq and 
increases strongly with halo mass and toward halo center. We find no detectable 
quenching enhancement for galaxies beyond ~ 2 i?vir around massive clusters once 
these galaxies have been decomposed into centrals and satellites. These trends imply 
that (1) galaxies experience no significant environmental effects until they cross within 
~ i?vir of a more massive host halo, (2) after this, star formation in active satellites 
continues to evolve in the same manner as active central galaxies for several Gyrs, 
and (3) once begun, satellite star formation quenching occurs rapidly. These results 
place strong constraints on satellite-specific quenching mechanisms, as we will discuss 
further in companion papers. 

Key words: methods: statistical - galaxies: groups: general - galaxies: clusters: 
general ~ galaxies: haloes - galaxies: evolution - galaxies: star formation 



1 INTRODUCTION 

Galaxy properties depend on their local environment. 
Galaxies in denser regions, such as groups and clusters, ex- 
hibit attenuated star formation rates (SFR), and therefore 
significantly higher red sequence fractions, as well as more el- 
liptical morphologies as compared with galaxies in less dense 
environments (Oemler||1974 Davis & Geller 1976 Dressier 



1980| [Postman fc Geller||1984| ). With the advent of voluml 
nous galaxy surveys such as the Sloan Digital Sky Survey 
(SDSS; York ot al. 2000| ), numerous works have quantified 
the correlations between these galaxies properties and their 
environment in considerable detail (e.g. 



vironmental impact on SFR (or color) is stronger than on 
morphology (Kauffma nn et ah] 2004 Blanton et al. 2005 
[Balfet al.„2008^ Bamford et al.J2009[ ): at fixed morphology, 
SFR still exhibits strong changes with environment, but at 
fixed SFR, there is almost no dependence of morphology 
on environment. Second, SFR/color and morphology pri- 
marily depend on small-scale (< 1 Mpc) environment, with 
little-to-no additional dependence on environment measured 
on larger scale s ([Hogg et al.||2004[ [Kaufi'mann et al.|r2004| 



Blanton et al. 20051. Several works have extended this re- 



sult through the use of galaxy group catalogs, showing that 
galaxy color and star formation history most directly depend 



Hogg et al. 



2004 



on the properties of their host dark matter halo (Blanton 



Balogh et aL][2004l [Kauffmann et al.||2Q04[ [Blanton et al.| [fc Berlind[[2007| [Wilman et al.[ [2010} [Tinker et ai.|[201l i' 



2005 1 , and more recent works have shown that similar envi- 

1 



ronmental dependence persists at least out to 



ciati et al.|2"006| [Cooper et al.|2007|[Peng et al.|2010 l 



(Cue- 



There are two particularly interesting features of the 
dependence of these properties on environment. First, en- 



Thus, understanding galaxy star formation requires identi- 
fying and understanding a galaxy's host halo properties. 

The idea that environmental dependence primarily 
manifests itself via host halo properties is motivated by the 
fact that the region within a halo is physically distinct from 
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the field environment. Tiie lialo virial radius corresponds to can merge witli one another (Makino & Hut 19971, a pro- 

a physical transition between infall and virialized motions, cess that may not be uncommon given that satellite orbits 

capable of supporting strong shock fro nts, at least in ha- are strongly correlated (Angulo et al. 2009 Wetzel et al. 

Ids m ore massive than ~ 10^^ M© (e.g., |Dekel fc Birnbo"im| |2009b|a[ [White et aL]|2010| |Cohn||2012[ ). These dynamical 



20061). Within the virial radius, high densities lead to strong 



tidal forces, and accreted gas is thermalized to high temper- 
ature and pressure. 'Satellite' galaxies that fall into more 
massive halos thus experience tidal and ram-pressure strip- 
ping, heating, and are unable to accrete gas as efficiently as 
in the field. Thus, within groups and clusters, it is important 
to distinguish satellite galaxies from the galaxy at the min- 
imum of the host halo potential well (the 'central' galaxy), 
which occupies a special dynamical location and is typically 
the oldest and most massive in the halo. 

In this paper, we focus on the SFRs of satellites galax- 
ies, how they depend on their host halo properties, and how 
they compare with central galaxies, which have not (yet) 
fallen into a more massive halo. Understanding satellite star 
formation is important not only for elucidating the physical 
processes that occur within groups and clusters, but also for 
a comprehensive understanding of galaxy evolution. For ex- 
ample, satellites are responsible for more than a 1/3 of the 
build-up of the red sequence population (van den Bosch et al. 



2008 Tinker fc Wetzel 2010). In addition, many methods of 



optically identifying galaxy groups/clusters rely on select- 
ing red sequence galaxies (e.g., Koester et al.|2007 l because 
they are more prominent within groups/clusters, making it 
important to characterize in detail how the satellite red se- 
quence fraction depends on halo mass and radius. 

Even though it has long been known that galaxies in 
groups/clusters are significantly more likely to have attenu- 
ated SFRs ([Dressier fc Gunn|1983||Balogh et al.|1997||Pog 



gianti et al. 1999), the mechanism(s) that dominates this 
transformation remains unclear, though several have been 
proposed. As a satellite galaxy orbits through its host halo's 
hot, virialized gas, ram-pressure can strip cold gas directly 
from the disc ( jGunn fc Gott|1972|[Abadi et al.|1999| ), caus- 
ing a rapid (< 500 Myr) quenching of star formation. There 
are striking, direct examples of ram-pressure stripping oc- 
curring for satellites in nearby clusters (e.g., Chung et al 
2009), some of which are even near the virial radius. While 



this suggests that ram-pressure stripping might be responsi- 
ble for satellite quenching in massive clusters, even at large 
radii, it remains unclear whether ram-pressure can act ef- 
ficiently in lower mass groups, particularly if virial shock 
fronts are not easily supported. 

Other mechanisms affect gas accretion, both into the 
subhalo that surrounds a satellite and cooling onto the 
disc. Referred to as 'strangulation' or 'starvation', this more 
gradual process involves the lack of accretion of new hot 
gas, as well as loss of existing extended hot gas via heat- 
ing/evaporation, tidal stripping, or ram-pressure stripping 



( Larson et al.|1980[|Balogh et al.|2000[[Kawata fc Mulchaey 



2008 McCarthy et al. 2008), causing SFR in the disc to 



decline on a cold gas consumption timescale. X-ray obser- 
vations show direct evidence for the truncation of hot gas 



halos around satellite galaxies in groups and clusters (Sun 
et al.|2007||Jehema et al.|2008[ ). 

Finally, galaxies in groups and clusters can interact with 
each other gravitationally. 'Harassment' involves high-speed. 



nearby encounters that tidally heat satellites ( jFarouki fc 
Shapiro|[l98l] [Moore et al.||1998[ ). More violently, satellites 



processes potentially could induce rapid gas consumption 
and thus quenching, in addition to driving morphological 
transformation. 

To understand which of the above physical process(es) 
dominates, it is crucial to have a detailed understanding of 
how satellite SFR depends on its host halo's mass and ra- 
dius within the halo, including whether there is a minimum 
halo mass for affecting satellite SFR. The above physical 
processes have different predictions for these dependences. 
However, the extent to which a satellite's SFR depends on its 
host halo mass, and how this compares with intrinsic depen- 
dence on its stellar mass, remains in debate. Several recent 
works using SDSS group/cluster catalogs found evidence 
that satellites are more likely to be quenched/red in more 
massive halos: Weinmann et al. (2006) examining both SFR 
and color using the group catalog finder of Yang et al. | ( [2005 ) ; 
Blanton fc Berlind[ (,2007) examining color using their own 



group catalog; Kimm et al. ([2009 1 examining UV-based SFR 



using the group catalog of [Yang et al. (20071; and Hansen 
et al. ( 2009 1 examining color in the MaxBCG catalog as 



a function of cluster richness. However, other works con- 
clude that the satellite SFR/color dependence on halo mass 
is non-existent or at least much smaller than the intrinsic de- 



pendence on stellar mass: Finn et al. (20081 examining SFR 



in the C4 cluster catalog as a function of cluster velocity 
dispersion; and van den Bosch et al. (2008)b and Pasquali 



et al. ( 2009 ) examining color and SFR, respectively, using 



the group catalog of Yang et al. (20071. Similarly, many re- 



cent works found that satellites are more likely to be red ( De 



Propris et aLl[2004l [Blanton fc Berlind 2007) [Hansen et al 



2009[) and have loweTSFRs (|Balogh et al.[[2000[ [EUingson 
et al.|2001|[Weinmann et al.|2006|[von der Linden et al.poiof 



Prescott et al.|20TT ) toward group/cluster center, including 
weak trends for groups at 2 ~ 1 (Gerke et al.|200 7') . Though 
again, there is debate about the steepness of radial gradi- 
ents, the extent to which they depend on halo mass, and 
whether they are a byproduct of satellite mass segregation 
(van den Bosch et al. 2008; Pasquah et al.[[2009l. 



There are several possible reasons for these disagree- 
ments. First, one must compare galaxies at fixed stellar mass 



because SFR and color depend strongly on mass ( Kauffmann 



et al. 2003), and more massive galaxies are more common 



in denser environments ( |Hogg et al. 2003) . The environmen- 
tal trends noted in many previous works were largely a re- 
flection of these intrinsic stellar mass dependencies. Second, 
some SFR metrics are more susceptible to systematic biases 
than others, and in particular, simple color cuts can over- 
estimate the quenched fraction because of dust reddening 
(e.g.. Mailer et al.|[2009 l. Third, in order to elucidate these 
trends, it is important to study satellite properties over a 
wide range of both satellite and host halo masses. 

Finally, to fully understand the physical mechanisms 
operating on satellites, it is crucial to go beyond simple 
SFR/color cuts and understand the nature of the full SFR 
distribution. Balogh et al. ( 2004 ) examined the u — r color 



distribution of galaxies in SDSS as a function of projected 
galaxy environment (defined by fifth nearest neighbor) and 
found that the color bimodality extended to all environ- 
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ments, and that the colors of blue galaxies do not depend 
significantly on environment. Using several different SDSS 



group catalogs, Skibba ( 2009 1 showed that the g~r color dis- 
tribution of satellites persists across group varying richness. 



More recently, McGee et al. (2011 1 found similar trends ex- 



amining the SFR distributions of all galaxies in groups and 
in the field both in SDSS and at z « 0.4. 

In this work, we explore the satellite galaxy SFR distri- 
bution in detail, including its dependence on host halo mass 
and halo-centric radius, to robustly characterize the prop- 
erties of satellite star formation evolution and quenching in 
groups and clusters. To observationally determine the host 
halo masses of SDSS galaxies, we construct galaxy group 



catalogs using a variant of the Yang et al. ( 2005 \ group- 



finding algorithm, which also allows us to examine how the 
SFRs of satellites differ from central galaxies of similar mass. 
Our galaxy sample and group finder allows us to explore a 
considerable range of halo and satellite masses to determine 
the regimes in which halo-specific processes are important. 
We examine galaxies selected in narrow bins of stellar mass 
to isolate host halo dependence from intrinsic stellar mass 
dependence, and we use spectroscopically-derived SFR mea- 
surements which are free from dust reddening effects. 

This paper represents the second in a series of four. In 
Tinker et al. ( 2011 1, hereafter referred to as Paper I we pre- 



sented our method for identifying galaxy groups in SDSS. 
Using the 400QA break, D„40Q0, as a measure of recent 
star formation history, we showed that central and satellite 
galaxy star formation is nearly independent of large-scale 
environmental density beyond its dark matter host halo. In 
this work, we explore the full SFR distribution of satellite 
galaxieajand its dependence on stellar mass, halo mass, and 
halo-centric radius to fully characterize the local density de- 



pendence of satellite star formation. In Wetzel et al. ( 2012a I , 



hereafter referred to as |Paper III| we will combine the results 
of this paper with a high-resolution cosmological simulation 
that tracks satellite infall times in order to examine satellite 
star formation histories and quenching timescales in detail. 



Finally, in Wetzel et al. (2012b I, hereafter referred to as Pa- 



|per IV| with the results of this paper we will use satellite 
orbital histories from simulation to test and constrain the 
physical mechanisms of satellite-specific quenching. 

For all calculations we use a fiat, ACDM cosmology of 
fini = 0.27, fib = 0.045, h = 0.7, n^ = 0.95 and ag = 
0.82, consistent with a wide array of observations (see, e.g., 
[Komatsu et al.|[2011[ and references therein). 



2 METHODS 

The details of our galaxy sample and group finding algo- 
rithm are discussed extensively in [Paper l| Here, we high- 
light the aspects relevant to this work. 



^ We have examined all trends in this paper selecting galaxies on 
both Mr and Mgtar as well as examining both D„4000 and SSFR, 
and while these lead to slight quantitative differences, they do not 
change any of our results qualitatively. 



2.1 Galaxy catalog 

To construct our galaxy sample, we use the New York Uni- 



versity Value-Added Galaxy Catalog ( Blanton et al.| 2005 I 
based on SDSS Data Release 7 ( Abazajian et al.|2009p . This 
spectroscopic catalog ensures that accurate galaxy redshifts, 
stellar masses, and SFRs can be obtained. Galaxy stellar 
masses are based on the kcorrect code of lBlanton fc RoweisI 



(20071, assuming a Chabricr (2003) initial mass function. 



To construct stellar mass-limited samples, we first con- 
struct two volume-limited samples of all galaxies with Mr < 
— 18 and Mr < —19, which contain galaxies out to 2; = 0.04 
and z — 0.06. Within these magnitude limits, we identify 
minimum complete stellar mass limits of 5 x 10^ M© and 
1.3 X 10^" M0, respectively. Combining these samples leads 
to an overall median redshift oi z = 0.045. 

For a galaxy star formation metric, we use specific star 
formation rate (SSFR), where SSFR = SFR/Af^tar. These 
values are based on the current releasaj of the spectral re- 

( 2004[ ), with updated pre- 



ductions of Brinchmann et al 



scriptions for active galactic nuclei (AGN) contamination 
and fiber aperture corrections following [Salim et al.| ( |2b07[ ). 
These SSFRs are primarily derived from emission lines 
(mostly Ha), but in the cases of strong AGN contamination 
or no measurable emission lines, the SSFRs are inferred from 
D„4000. Roughly, SSFR > 10"" yr"^ are based almost en- 
tirely on Ha, 10"^^ < SSFR < 10"" yr"^ are based on a 
combination of emission lines, and SSFR < 10~^^ yr~^ are 
based almost entirely on D„4000 and should be considered 
upper limits to the true value. 

These SSFRs have also been corrected for aperture bias 
because the SDSS spectral fiber size of 3 arcsec (1.5 h~^ kpc 
at our median redshift) can be smaller than a galaxy's ef- 
fective radius. Specifically, [Brinchmann et al.| ( |2004[ ) mea- 
sured the distribution of SSFRs at a given g — r and r — i 
color within the fiber for all galaxies and applied this empir- 
ical calibration to the photometry outside the fiber. Their 
updated, aperture-corrected SSFRs yield good agreement 
with full-galaxy aperture, UV-hased SSFR measurements 
of |Salimet ar] ( [2057t . 

Galaxy pairs closer than 55 arcsec (30/i~^kpc at our 
median redshift) are too close to both receive spectral fibers 
('fiber collisions'), so 5% of the galaxies in our sample do not 
have spectra. These galaxies are assigned the same redshift 
as their nearest non-collided galaxy, and we assign SSFR by 
sampling randomly from the distributions of non-collided 
galajcies with the same r-band magnitude and g — r color 
within the volume-limited sample. We have analyzed all of 
our results both using and excluding fiber-collided galaxies, 
and both produce consistent results, though to be safe in 
examining trends with radius we only examine radial bins 
where the fiber-collided fraction is less than 10%. 

We emphasize that the use of spectroscopically-derived 
SSFRs is critical for our analysis because dust reddening can 
severely overestimate the quenched population when using 
red/blue color cuts, particularly at lower mass where gas 
and dust fractions are higher. For our lowest mass galaxies 
(Mstar < 10^° M0) 37% of red {g-r> 0.76) galaxies are in 

, and 21% of active galaxies 



fact active (SSFR > 10" 

appear red (see Fig. 2 in Paper II. We have confirmed that 



yr 



http : //www . mpa-garching .mpg . de/SDSS/DR7/ 
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this population of active but red galaxies is caused primar- 
ily by dusty spirals by applying Galaxy Zoo morphological 
classifications (Lintott et al. 2011 1 to this population: ~ 70% 
are identified as spirals, and ~ 50% are edge-on spirals (see 
also [Masters et"al]|2010 1 . 



2.2 Galaxy group catalog 

Motivated by the paradigm that all galaxies reside in host 
dark matter halos, we identify the halo properties of galax- 
ies by using a modified implementation of the group-finding 
algorithm detailed in lYang et al.l (120051 120071). We define 



a galaxy group as a set of galaxies that occupy the same 
host dark matter halo, and we will use 'group' and 'halo' 
interchangeably henceforth. 

For our group catalog, we define dark matter halos 
such that the mean matter density interior to the virial 
radius is 200 times the mean background matter density: 
Af200m = 200pm^7'"-R200m - Usiug the host halo mass func- 
tion from Tinker et al. ( 2008 1 combined with the satel- 



lite subhalo mass function from Tinker & Wetzel (20101, 



we first assign a tentative (sub)halo mass to each galaxy 
by matching the abundance of objects above a given mass: 
n(> M(sub)haio) = n{> Mstar). Each galaxy then has an 
associated (sub)halo virial radius and velocity dispersion. 
Starting with the most massive galaxy in the sample, we 
assign to each nearby, lower mass galaxy a probability of 
being a satellite member of its group using a matched filter 
in both radius and line-of-sight velocity difference. We re- 
peat this process for less massive galaxies on down the mass 
function, skipping those that have already been assigned as 
satellites in a group, until each galaxy is assigned to a (tenta- 
tive) group. Using now the total stellar mass of each group, 
we update the halo mass of each group based on abundance 
matching, but now using only the host halo mass function 
(ignoring subhalos): n(> Mhaio) = n(> Mstar.group). With 
these updated halo masses, we re-assign satellite member- 
ship probabilities, we update the halo masses again, and 
we iterate until convergence. Every group thus contains one 
'central' galaxy, which by definition is the most massive, and 
a group can contain zero, one, or many satellites. 

We construct groups at z < 0.04 using the Mstar > 
10^'^ Mq sample and at 2: = 0.04 - 0.06 using the Mstar > 
10^°-^ Mq sample. This leads to ~ 19000, 2200, 160 groups 
with M2oom > 10", 10^^, 10^" M0, respectively, up to the 
most massive group at 10^^ M0 hosting 269 satellites. Our 
catalog allows us to examine satellites with good statistics 
in the mass range Mstar = 5 x l(f — 2 x 10^^ Mq. 

We define a group's center by the location of its cen- 
tral (most massive) galaxy. However, in reality the most 
massive galaxy is not always at the minimum of a host 
halo's potential well, for two reasons. First, the galaxy at 
the minimum of the potential might not be the most mas- 



sive in the halo. Using the Yang et al. (20071 group catalog, 
[Skibba et al.| ( |2011[ ) found evidence that up to 40% of mas- 
sive (> 10 Mq) halos have galaxies closest to halo cen- 
ter that are not the brightest, with a decreasing fraction at 
lower halo mass. However, these are typically cases where 
two galaxies in a halo have similarly high luminosities, and 
in this regime galaxies exhibit highly quenched SSFRs al- 
most regardless of central vs. satellite demarcation, so we 
do not expect that this effect will bias our satellite SSFR re- 



sults significantly. Second, the (true) central galaxy can be 
offset from the minimum of the halo potential as a result of 
dynamics processes such as merging. For example, for rich 
clusters in the SDSS MaxBCG catalog with X-ray detection, 
the median offset between the brightest cluste r galaxy and 
X-ray center is 58 /i~^ kpc ( Koester et al.|2007 l. To test the 
impact of both of these possible biases, we have examined 
our radius trends centering on both the most massive galaxy 
and the luminosity- weighted center of the group, and we find 
no significant difference in our results. 

Because of projection effects and redshift-space distor- 
tions, some central galaxies are inevitably misassigned as 
'satellites' of higher mass halos (reducing group purity), 
and conversely some satellites are misassigned as 'central' 
galaxies of lower mass halos (reducing group completeness) . 
As detailed in [Paper l| we find that an average of ~ 10% 
of galaxies have been misassigned, with a stronger effect 
for more massive satellites in lower mass halos (see also 
Yang et al.||2007l. Because central galaxies exhibit a lower 



quenched fraction than satellites, and because they outnum- 
ber satellites, the primary resultant biasing is an underesti- 
mation of the quenched fraction for satellites. This effect also 
causes an overestimation of the quenched fraction for central 
galaxies, but because central galaxies outnumber satellites 
this effect is relatively small. In this work, we do not attempt 
to correct our results for these biases given their non-trivial 
dependence on halo mass and radius. 

Finally, we have checked our results against a SDSS 
DR7 version of the |Yang et al.| ( [2007[ ) group catalog (kindly 
provided by Frank van den Bosch), incorporating the same 
kcorrect stellar masses and SSFR measurements we use 
here but retaining differences in assumed cosmology, halo 
virial definitions, and group-making methodology. All of the 
results are consistent within errors. 



3 MASS DEPENDENCE 

3.1 Star formation rate distribution 

We begin by exploring the galaxy SSFR distribution and 
its dependence on galaxy and host halo mass. Fig. [l] shows 
the distribution of SSFRs for galaxies decomposed into cen- 
trals and satellites. Each panel shows a different stellar mass 
bin, and satellites are divided into halo mass bins. Examin- 
ing first the central galaxy distribution, at all stellar masses 
it exhibits a clear bimodality with a break at 10~^^ yr"^- 
More massive central galaxies have a decreased likelihood of 
high SSFRs and an enhanced likelihood of low SSFRs, but 
the break location remains fixed. This SSFR bimodality and 
fixed break value have been noted by several authors (e.g., 
Brinchmann et al.|2004 Kauffmann et al.|2004) . 

Satellites exhibit a clear deficit/excess at high/low 
SSFR as compared with central galaxies of the same stellar 
mass. This difference is smaller in lower mass halos, but we 
emphasize that Fig. Ilk shows that a clear difference persists 
even at Mhaio J; 10 Mq, and that this is a lower limit to 
the actual difference because we have not corrected for pu- 
rity and completeness effects. Thus, we find no evidence for 
a minimum halo mass for affecting satellite star formation. 
In more massive halos, the relative fraction of high SSFR 
galaxies decreases monotonically with halo mass across all 
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Figure 1. Specific star formation rate (SSFR) distribution for 
central galaxies (solid) and satellites in bins of host halo mass 
(dashed and dot-dashed) . All galaxies in all halo masses we probe 
exhibit similar bimodal SSFR distributions. Satellites in more 
massive halos are more likely to be quenched, but the SSFR values 
of the high SSFR peak, the break at SSFR fs 10"" yr"\ and 
the fraction of galaxies at the break remain unchanged. Dotted 
curves in (a) show example fits to a skewed double Gaussian. 
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Figure 2. Active galaxy SSFR distribution peak vs. host halo 
mass, in bins of galaxy stellar mass, as determined by fitting 
a skewed, double Gaussian to each full SSFR distribution (see 
Fig.llVi for example). Points show central galaxies, plotted at the 
median host halo mass of each sample, and curves show satellites. 
While more massive galaxies show a trend toward lower SSFR 
peak values, there is no significant dependence on host halo mass. 



halo masses we probe. Remarkably, though, the SSFR values 
where the satelUte bimodaUty break and active peak occur 
do not change, nor does the fraction of satellites lying at the 
10~^^ yr~^ bimodality break (the 'green valley'). 

We examine the invariance of the high SSFR peak more 
quantitatively by fitting a skewed, double Gaussian to the 
full SSFR distributions in Fig. [T] with example fits shown 
by dotted curves. Fig. [2] shows the peak of the distribution 
for high SSFR ('active') galaxies as a function of host halo 
mass. (Note that SSFR = lO^^'''^ yi'~^ corresponds to a 
rate that doubles mass in a Hubble time.) Square points 
at left show central galaxies, plotted at the median halo 
mass corresponding to the stellar mass bin (recall that cen- 
tral galaxies have a tight relation between stellar and halo 
mass), while curves show satellites. More massive galaxies 
exhibit systematically lower SSFR peak values, but within 
a given stellar mass bin there is no significant change with 
halo mass. Thus, the SFRs of galaxies that remain active 
are not affected after falling into a host halo. We discuss the 
significance of these SSFR bimodahty results further in i|5] 

We do not examine the peak at low SSFR, as it is 
partially an artifact of the spectral reductions of |Brinch-| 
|mann et al.| ( |2004| ), where low SSFR galaxies with no de- 
tectable emission lines are assigned SSFRs based on D„4000, 
which effectively leads to an upper limit value of SSFR ~ 
10~^^ yr~^. Thus, while the rise of the distribution below 

^^ yr~^ is real, the strong sharpness of the peak near 
12 ' „-i 



10 

10^^^ yi'^^ is largely artificial, and in reality the distribu- 
tion exhibits a tail to much lower SSFR ( jSalim et al. 12007] ). 
Given the invariance of the value of the SSFR bimodal- 
ity break across galaxy mass, host halo mass, and satellite- 
central demarcation, we divide galaxies into those with 
SSFR > 10^^^ yr~^) referred to as 'active', and those with 
SSFR < 10"" yr"\ referred to as 'quenched'. With this 
active/quenched demarcation, we will examine the fraction 
of galaxies at a particular stellar mass that have quenched 
SSFRs, that is, the 'quenched fraction'. This quantity rep- 
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resents a more physically meaningful version of the 'red se- 
quence fraction', and the two are effectively synonymous in 
the limit of no dust reddening. 

3.2 Quenched fraction 

We now explore the galaxy and host halo mass dependence 
of the quenched fraction in detail. Fig.lsk shows that the like- 
lihood of central galaxies to be quenched increases strongly 
with stellar mass, rising from 20% to almost 95% across a 
factor of 15 in stellar mass. This mass range represents the 
strongest transition in galaxy star formation properties (e.g., 
Kauffmann et al.|2003 l. By contrast, the satellite quenched 
fractions rise more gradually with stellar mass at fixed halo 
mass. This is because low-mass satellites are more likely to 
be quenched, leaving less room for stellar mass dependence. 

Fig. [3b shows the same but demonstrates more directly 
how the quenched fraction increases with halo mass at fixed 
stellar mass. Square points show central galaxies plotted at 
their median host halo mass for each stellar mass bin, while 
curves show satellites. Satellites reside in higher mass host 
halos than central galaxies of the same stellar mass because, 
by definition, satellites are not the most massive galaxy in 
their group, so the total stellar mass, and hence halo mass, of 
a satellite's group is typically higher than if it were a (most 
massive) central galaxy. At a given stellar mass, satellites 
are always more likely to be quenched than central galaxies. 
This is true even in the lowest halo mass scales we probe 
(3 X 10^^ Mq), again implying no minimum halo mass for 
satellite quenching. Comparing instead satellites to central 
galaxies in the same halo masses, central galaxies are always 
more likely to be quenched than the satellites in their same 
host halo. This latter trend arises from the strong depen- 
dence of the quenched fraction on stellar mass. 

Fig. [3b shows that, at fixed stellar mass, transition- 
ing from a central to a satellite galaxy in a progressively 
more massive host halo results in a significant increase in 
likelihood of quenching, but this increase with halo mass is 
weaker for more massive galaxies (see Pasquali et al.[[2010 
for a similar trend for galaxy age). At face value, this might 
seem to imply that more massive galaxies are less affected by 
satellite-specific processes. However, one must take care in 
interpreting absolute changes in quenched fraction. At lower 
stellar mass, central galaxies are less likely to be quenched 
and thus are able to experience a larger absolute change in 
quenched fraction after infall. At higher stellar mass, there 
is less room for additional satellite-specific quenching. 

We can examine host halo mass dependence more di- 
rectly by removing the intrinsic dependence of the quenched 
fraction on stellar mass. Specifically, we examine the relative 
change in quenched fraction between central and satellite 
galaxies of the same stellar mass, and we scale this by the 
fraction of active (able to be quenched) central galaxies: 
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where fg and /a are the fraction of quenched and active 
galaxies, respectively, related via /q = 1 — /a and measured 
independently for central and satellite galaxies. A similar 



statistic was used in van den Bosch et al. ( 2008 1 and Tinker 



fc Wetzel (20101, and more recently by Peng et al. (20111, 
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Figure 3. (a): Galaxy quenched fraction (SSFR < 10 



yr 



-^) 



vs. stellar mass, for central galaxies (solid) and satellites in bins 
of host halo mass (dashed and dot-dashed). Both central and 
satellite quenched fractions increase strongly with stellar mass. 
(b): Same, but as a function of host halo mass split into bins 
of stellar mass. Curves show satellites while points show central 
galaxies at their median host halo mass. The satellite quenched 
fraction increases with halo mass and is always higher than the 
central galaxy quenched fraction at the same stellar mass, (c): 
Satellite quenched fraction excess, given by equation ([Til, vs. host 
halo mass. This excess is nearly independent of satellite mass. 
Error bars show 68% confidence interval for a beta distribution 
(see |Cameron|201 1[| , and curves are offset slightly for clarity. 



but our interpretation of it differs from those. If the central 
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galaxy active fraction has not evolved since the time of satel- 
lite infall, such that /^"(zinf) = /^"(^ « 0), then /5^Lce== 
represents the fraction of satellites that quenched after infall. 
However, as we will show in [Paper III| the typical redshift 
of infall for current satellites in our mass range is z « 0.5, 
so the assumption of no central galaxy evolution since the 



time infall is not a good one (see also Tinker & Wetzel 2010 1 



To the extent that the central galaxy quenched fraction has 
increased over time, /q'cxccss more accurately represents the 
excess fraction of satellites that quenched after infall that 
would not have quenched had they remained central galax- 
ies, as can be seen by the right-hand side of equation ([l]). 

Fig.™ shows that /q'excess depends only weakly on stel- 
lar mass at low halo mass, and it is remarkably independent 
of stellar mass in halos > 3 x 10^^ M©. (While the error at 
high stellar mass is large, the trend is consistent with that of 
lower stellar mass.) Moreover, the dependence of /g'excess on 
stellar mass in lower mass halos does not necessarily imply 
a change in satellite quenching efficiency, for two reasons. 
First, satellite dynamical infall times are shortest when the 
satellite-to-halo mass ratio is closest, being shorter than a 
Hubble time when the satellite halo to host halo mass ra- 
tio is closer than ~ 0.1 ( Boylan-Kolchin et al.pOOS Jiang 



et al. 2008 Wetzel & White 20101. Thus, surviving satel 



lites in this regime have necessarily fallen in more recently 
and have not had as much time to quench. Second, purity 
and completeness in the group catalog are worse when the 
satellite-to- halo mass ratio is high, so satellite quenched frac- 
tions are most underestimated in this regime. 

These satellite quenched fraction excess results show 
that the apparent reduction of the dependence of the 
quenched fraction on host halo mass at higher stellar mass 
primarily reflects that a higher fraction of more massive 
satellites quenched prior to infall. 



4 RADIAL DEPENDENCE 

4.1 Quenched fraction within R^ir 

Having explored the mass dependencies of the SSFR distri- 
bution and quenched fraction, we next examine how these 
quantities depend on location within the host halo. Exam- 
ining the dependence on halo-centric radius is particularly 
informative: because dynamical friction causes a correlation 



between satellite radius and time since infall iGao et al. 
|2004,_Weinmann et al.|2011[ also |Paper IVp , radial gradients 
translate into an evolutionary sequence of star formation. 

In examining satellite radial gradients, one must ensure 
that the results are not affected by possible mass segrega- 
tion (more massive satellite preferentially being at smaller 
radii). We find no evidence for satellite mass segregation 
at any halo mass in our group catalog, in agreement with 
previous works ([Pracy et al.||2005| [Hudson et al.||2010[ |von 



der Linden et al. 20101, but see also van den Bosch et al 



(|2008)b. Furthermore, we examine satellites in narrow bins 
of stellar mass, thus minimizing the effects of any possible 
mass segregation. 

Fig. |4k and b show the satellite quenched fraction as 
a function of projected radius scaled to the host halo virial 
radius, varying both halos mass and satellite mass, respec- 
tively. For comparison, points show the average quenched 
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Figure 4. Galaxy quenched fraction vs. projected radius, scaled 
to the host halo virial radius. Curves show satellites and points 
show central galaxies, regardless of location, of the same stel- 
lar mass, (a) and (b): Varying halo mass and satellite mass, 
respectively. The quenched fraction increases with decreasing ra- 
dius, and at all radii the satellite quenched fraction is higher in 
more massive halos and for more massive galaxies, (c): Satellite 
quenched fraction excess, given by equation (IlJ, is independent of 
satellite mass at all radii. Curves are offset slightly for clarity, and 
the highest stellar mass has reduced radial binning given limited 
statistics. 
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fraction for all central galaxies of the same stellar mass. At 
all satellite and halo masses, the quenched fraction increases 
toward smaller radius. These gradients are strong, particu- 
larly for lower mass galaxies, which exhibit a 4 times higher 
quenched fraction at 0.05 -Rvir than in the field. These gra- 
dients imply that a satellite's quenching likelihood depends 
on its time since infall and/or whether it passed through 
denser regions of its host halo. 

Interestingly, the quenched fraction is higher in more 
massive halos even at fixed R/Rvir, by an amount that 
is independent of radius. This result implies that satellite 
quenching does not simply depend on its currently observed 
local density, because at a given R/Rvii, halos of any mass 
have the same density (modulo minor dependence on halo 
concentration). This does not rule out the possibility that 
satellite quenching could occur via crossing above some fixed 
background density threshold, but it means that orbital evo- 
lution as/after quenching occurs changes the observed corre- 
lation with density. In other words, currently observed local 
density is not necessarily indicative of the maximum local 
density a galaxy has ever experienced. Thus, understand- 
ing satellite quenching requires understanding the detailed 
orbital histories of satellites, as we will examine in |Paper IV| 

The results of Fig. [4V) show that radial gradients are 
shallower at higher stellar mass, which again might seem 
to imply that satellite quenching efficiency decreases with 
satellite mass. To examine this more carefully, we again mea- 
sure the satellite quenched fraction excess, given by equation 
([I|, in bins of scaled radius. Fig. En shows that, remarkably, 
/q'cxcgss is independent of satellite mass at all radii. (Again, 
the error bars at high stellar mass are large, but the trend 
is consistent.) This result strongly supports our conclusion 
that the apparent reduced dependence on local density at 
higher galaxy mass simply reflects a higher fraction of satel- 
lites having quenched prior to infall. 

Note that the results in Fig. |4] are for projected radius 
and incorporate purity and completeness effects. At small 
projected radii the primary contamination is the incorpora- 
tion of satellites in the same halo but at larger (line-of-sight) 
radii, while at large radii the primary contamination is from 
other central galaxies scattering into the halo. Both effects 
lead to an underestimation of the true satellite quenched 
fraction at a given radius, so the true 3-D radial gradients 
could be somewhat steeper or shallower than in Fig. [I] 



4.2 Quenched fraction outside i?vir 

The results from Fig. |4] show that the quenched fraction 
for satellites can be significantly higher than that of central 
galaxies even for satellites near the virial radius in high- 
mass halos. This effect may be driven in part by these 
satellites having passed within much smaller radii, given 



that many satellites are on highly elliptical orbits (Ben- 
|son|2005|[Khochfar fc Burkert|2006} [Wctzcl 20fT|, and that 
galaxy groups falling into more massive halos can bring in 
already quenched satellites (Berrier et al.|[2009[ |Li et al 



2009 McGee et al.|[2009K However, it is interesting to con- 



out to 2 — 4 _Rvir around massive clusters ( Hansen et al.|2009 
von der Linden et al.|2010 l. 

Fig. [5] shows the galaxy quenched fraction as a func- 
tion of projected radius, similar to Fig. Ill except including 
all galaxies (centrals and satellites) within a projected ra- 
dius of 10 -Rvir around massive clusters, regardless of whether 
they are members of the selected halos or reside in different 
halos. (Because we impose no redshift cut here, some galax- 
ies within _Rvir are not members of the selected cluster-mass 
halos.) The quenched fraction for all galaxies (solid curve) 
exhibits a clear enhancement out to ~ 3 — 4_Rvir. However, 
if we decompose this population into central and satellite 
galaxies according to the group catalog, it is clear that the 
enhancement beyond Rvir is driven primarily by satellites in 
different halos and a declining ratio of satellites-to-centrals 
as a function of radius from the clusters. The latter sim- 
ply reflects an increased likelihood of finding massive groups 
closer to clusters. Indeed, compared with the central galaxy 
quenched fraction measured at 5 — 10 7?vir (dotted curve), 
central galaxies deviate significantly from the field value only 
within ~ 2_Rvir- At these radii, the enhancement is likely 
caused by (1) massive halos being highly ellipsoidal, and (2) 
satellites on orbits with apocenters beyond the virial radius 
( |Gill et al.||2005| ILudloweFal. 2009 Wang et aH2009|). We 

jPaper IV| 



explore the latter explanation in more detail in 



See also Y. Wang et al. ( 2009 1 for similar trends for much 



sider whether environment affects star formation in galax- 
ies beyond the virial radius by examining whether the ra- 
dial gradients of Fig. [4] persist to larger radii. An enhanced 
red/quenched fraction has been noted in examining galaxies 



fainter galaxies (M^ < —17). 

Thus, by decomposing galaxies beyond R^ii into cen- 
trals and satellites, we find no compelling evidence that cen- 
tral galaxy star formation is affected significantly prior to 
infall. Instead, the SSFR gradients beyond the virial radius 
are caused primarily by satellites in groups and that massive 
groups are more likely to be found closer to clusters. This 
agrees with our results in |Paper"T| where we showed that the 
change in the overall galaxy quenched fraction with large- 
scale (10/i~^Mpc) environmental density is caused simply 
by a changing halo mass function. Thus, all environmental 
effects on galaxy star formation are consistent with being 
governed only by a galaxy's host halo. 



4.3 Star formation rate distribution 

Finally, to understand more fully the radial gradients in 
Fig. [4] we examine how the full SSFR distribution changes 
with halo-centric radius. Fig. |6] shows the satellite SSFR 
distribution in bins of scaled radius, with panels showing 
different host halo masses. (We show low-mass satellites, for 
which changes in the SSFR distribution are most clearly vis- 
ible, though similar trends persist for higher mass satellites.) 
For comparison, the solid curves show central galaxies of the 
same stellar mass. Remarkably, the bimodal nature of the 
SSFR distribution persists at all radii within all halo masses 
we probe. Similar to the halo mass trends in Fig. [21 there 
is no significant change in the active galaxy SSFR peak lo- 
cation with radius at the ~ 0.1 dex level, nor is there a 
significant change in the location of the bimodality break or 
the fraction of galaxies there. Thus, star formation in cur- 
rently active satellites is not affected by local density or time 
since infall. We now discuss the strong implications of this 
persistent SSFR bimodality. 
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Figure 5. Galaxy quenched fraction vs. projected radius, scaled to the host halo virial radius, including all galaxies within 10-R200m 
(projected) regardless of group assignment. Panels show different stellar mass bins around cluster-mass halos (M200m > 10^'* Mq). Solid 
curves show all galaxies, points show galaxies split into satellites and centrals, and dotted curve shows the average value for central 
galaxies at 5 — 10 -R200m- Central galaxies deviate significantly from the field value only within ~ 2 -R200m- 



5 IMPLICATIONS OF THE STAR 
FORMATION RATE BIMODALITY 

We have shown that all galaxies, central and satellite, ex- 
hibit a similar bimodal distribution of SSFRs, with a deficit 
at SSFR ^ 10~^^ yr^^ and the same active galaxy SSFR 
peak at fixed stellar mass. This bimodality extends across all 
stellar masses, host halo masses, and halo-centric radii; in no 
regime is there ever a pile-up of galaxies at intermediate SS- 
FRs in the 'green valley'. Thus, halo dependence of satellite 
galaxy SFR manifests itself simply as a decrease/increase 
in the fraction of active/quenched satellites with increasing 
halo mass and with decreasing halo-centric radius. 

This 'universal' SSFR bimodality places strong con- 
straints on satellite-specific quenching mechanisms and 
timescales. In particular, to the extent that halo-centric ra- 
dius selection represents an evolutionary sequence, the re- 
sults of Figs. IT] and [6] have three clear implications. First, 
the fact that active satellites have the same SSFR peak and 
distribution as active central galaxies, even within 0.17?vir, 
implies that a significant fraction of satellites have evolved in 
the same manner as central galaxies, unaffected by host halo 
processes after infall. This suggests that any halo-specific 
quenching process takes significant time (at least an orbital 
time of several Gyrs) to start to affect satellite star forma- 
tion (see also McGee et al.|2011 1. 

Second, the persistent bimodality break implies that 
once quenching begins, the transition to quenched SSFRs 
always occurs on a rapid timescale. Specifically, once be- 
gun, the fading timescale must be significantly shorter than 
the timescale over which satellites remain active after in- 
fall (several Gyrs), and likely shorter than molecular gas 
consumption timescale (~ 2.4 Gyr, e.g., Bigiel et al.||2011l, 



timescale after infall, we would observe an excess of galaxies 
at intermediate SSFR values in the 'green valley'. However, 
the fact that the fraction of satellites at SSFR « 10~^^ yr~^ 
is both non-zero and constant implies that satellite SFR fad- 
ing does take cosmologically non-trivial time to occur, and 
that this fading timescale does not depend on host halo mass 
or halo-centric radius. It remains unclear if this scenario can 
be achieved merely by a sharp cold gas density threshold for 



sustaining star formation (see e.g., Krumholz et al. 20091 



or if some additional process is needed to accelerate cold gas 
consumption and SFR decay. We will return to this question 
in more detail in [Paper III[ 

Third, the deviations of the satellite SSFR distribution 
from that of central galaxies of the same stellar mass down 
to the lowest halo mass scales we probe (3 x 10^^ Mq) means 
that we find no minimum halo mass for satellite-specific pro- 
cesses. This result is consistent with ToUerud et al. (20111, 



who found an enhanced red fraction of Large Magellanic 
Cloud-mass galaxies if they are near Milky Way-mass galax- 
ies in SDSS. Furthermore, the increasing quenched fraction 
with halo mass implies either that satellite quenching is more 
efficient and begins more quickly in higher mass halos, or 
simply that satellites in more massive halos have preferen- 
tially been satellites for longer, possibly from satellite pre- 
processing in groups. In [Paper lH] and [Paper IV} we will 
demonstrate that the latter explanation is more likely. 

Thus, the persistent bimodality implies that a delayed- 
then-rapid quenching mechanism must drive satellite 
quenching. While the implied long timescale for satellite 
quenching has led several authors to propose strangulation 



as the dominant mechanism (e.g., Balogh et al. 2000 McGee 



because if satellite SFR simply were to fade slowly over this 



et al.|2009 Weinmann et al.|2010 l, it is unclear that a mode 
of strangulation in which star formation fades slowly upon 
infall can retain the SSFR bimodality at all radii without 
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Figure 6. SSFR distribution for central galaxies (solid) and 
satellites in bins of projected radius, scaled to the host halo virial 
radius (dashed and dot-dashed). Panels show different host halo 
mass bins. The SSFR, bimodality persists at all radii, with a con- 
stant location for the active galaxy peak, bimodality break, and 
fraction of galaxies at the break. The only difference is a decreased 
fraction of active galaxies at smaller radii. 



putting too many galaxies at intermediate SSFR (though see 
Weinmann et al.|[20To| and also ( [Balogh et aL|[2009l )). One 
possibility is that, while the satellite's hot gas halo is being 
stripped from the outside, the innermost hot gas remains 
unaffected such that cooling and star formation continue as 
normal for several Gyrs. This idea is consistent with the 



X-ray observations of Sun et al. (20071 and Jeltema et al 



( 2008 1 , who found that more than half of bright satellites 



have significantly truncated yet detectable hot gas halos. 

A natural mechanism to rapidly quench star formation 
on a short timescale (< 500 Myr) is ram-pressure strip- 
ping, which can remove cold gas directly from a satellite 
galaxy's disc. However, while ram-pressure is observed to 
act on galaxies in high-mass clusters, it has not been ob- 
served in lower mass groups, where lower halo gas densities 
and satellite velocities likely lower its efficiency. Moreover 
halos of Mvir J; 10^^ M© are not expected to have virial 
shock fronts which support hot, virialized gas within the 
halo ( Dekel fc Birnboim||2006 l, so in this mass regime it is 
not clear that either strangulation or ram-pressure can be 
efficient. This may suggest the need for tidal stripping, or 



for harassment and/or mergers induce rapid cold gas con- 
sumption that quenches star formation on short timescales. 
In |Paper IV[ we will test quantitatively whether these mech- 
anisms can reproduce the galaxy mass, halo mass, and halo- 
centric radius dependencies that we have outlined here. 



6 SUMMARY AND CONCLUSION 

Using galaxy group catalogs created from SDSS Data Re- 
lease 7, we examined in detail the SSFRs of satellite galaxies 
and how they depend on galaxy mass, host halo mass, and 
halo-centric radius. Our galaxy sample and group catalog 
provide good statistics across of wide range of masses both 
for satellites (Mstar = 5 x 10^ - 2 x lO" M©) and their 
host halos (M2oom = 3 x lO" - 10^^ M©), and our use of 
spectroscopically determined SSFR means our results are in- 
sensitive to dust reddening, important for low-mass galaxies 
Our main results are as follows. 

Persistent Bimodality: All galaxies that we probe, 
regardless of galaxy mass, host halo mass, and halo-centric 
radius, exhibit similar bimodal SSFR distributions, includ- 
ing a fixed break at SSFR = lO"'^^ yr~^, the same fraction 
of galaxies at the break, and the same SSFR peak for ac- 
tive galaxies. Satellites are simply more likely to lie on the 
quenched side of the distribution. This persistent bimodality 
implies that a significant fraction of satellites have evolved 
in the same way as central galaxies, independent of their 
host halo, for several Gyrs. Once begun, the satellite SFR 
fading process is rapid as compared with time over which 
satellites remain active after infall and thus is likely shorter 
than a ~ 2.4 Gyr cold gas consumption timescale, but it 
does not depend on host halo mass or halo-centric radius. 

Clear host halo mass and radius dependences 
within the virial radius: At fixed stellar mass, satellites 
are more likely to be quenched than central galaxies, even 
down to the lowest halo masses we probe (3 x 10^^ Mq). 
This quenched likelihood increases at smaller halo-centric 
radius, further implying that satellite quenching takes sev- 
eral Gyrs to occur after infall. The quenched likelihood also 
monotonically increases with halo mass across our entire 
halo mass range, even at fixed R/Rvh, which implies that 
satellite quenching does not simply depend on current local 
density. It is not immediately clear whether satellite quench- 
ing is more efficient/rapid in higher mass halos, or surviving 
satellites in that regime have simply been satellites longer. 

Unclear stellar mass dependence: The central 
galaxy quenched fraction increases significantly with stellar 
mass. While more massive satellites show a weaker absolute 
change in quenched fraction with host halo mass or radius, 
this is largely a manifestation of central galaxy mass depen- 
dence prior to infall. The satellite quenched fraction excess, 
which more directly measures quenching efficiency, shows 
little-to-no dependence on galaxy mass. 

Importance of satellite-specific evolution: In |Pa-| 
|per I[ we showed that the star formation histories and 
quenched fractions of central galaxies at fixed galaxy (or 
halo) mass are essentially independent of their large-scale 
(10/i~^Mpc) environment. Here, we have shown that this 
lack of environmental dependence extends to central galax- 
ies as close as ~ 2 Rvh to massive clusters. Together, these 
results demonstrate the importance of satellites in under- 
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standing galaxy evolution, because satellite-specific pro- 
cesses are the only significant environmental processes that 
affect galaxy star formation and color. 

The results of the satellite quenched fraction excess sug- 
gest that infalling active satellites quench at least as effi- 
ciently/rapidly at high stellar mass as at low stellar mass. 
However, it remains true that satellite-specific processes are 
less important for building up the red sequence at higher 
stellar mass, because more galaxies quench as central galax- 



ies prior to infall in that regime (see also van den Bosch 



et al. 2008 1. Unfortunately, it is difficult to determine di- 
rectly from fcflxccss whether satellite quenching is fully in- 
dependent of galaxy mass because doing so requires knowing 
in detail the quenched fraction as a function of stellar mass 
prior to satellite infall. In [Paper IIll we will use a high- 
resolution cosmological simulation to obtain the infall time 
distributions of satellites, combined with observational con- 
straints on high-redshift quenched fractions, to examine this 
issue more carefully. 

The dependences on mass and halo-centric radius that 
we see broadly agree with many previous works using SDSS 
catalogs. This includes a monotonic increase in quenched 
fraction with halo mass mass which nearly equals the in- 



trinsic dependence on galaxy mass iWeinma nn et al.||2006[ 
Blanton fc Berlind|[2007l [Hansen et al. 2009! [Kimm et al'[ 
20091, as well as satellites being more likely to be quenched 
at smaller radii (Weinmann et aLl|2006 Blanton & Berlind 



2007 Hansen et al. 2009 1 with less massive satellites exhibit- 



ing stronger radial gradients (von der Linden et al. 20101 



Our catalog enabled us to extend many of these results to a 
wider range of galaxy and halo masses, particularly at the 
low mass end, and our use of spectroscopically-derived SSFR 
ensures that of quenched fraction results are not biased by 
dust in that regime. Furthermore, our SSFR bimodality re- 
sults qualitatively agree with previous results on the color 



bimodality as a function of projected galaxy density ( Balogh 



et al. 2004 Baldry et al. 



McGeeet al.'(2011) and 



2006 1 and the more recent results of 



Peng et al. ( 2011 1, who showed sim- 



ilar SSFR distributions and active galaxy SSFR values when 
splitting all galaxies into groups vs. the field. We have shown 
that this 'universal' SSFR bimodality extends to specifically 
satellite galaxies in all halo masses and at all radii. 

Using color cuts to measure red galaxy fractions, [van[ 



den Bosch et al. ( 2008 1, and more recently Peng et al. ( 2011 1, 



employed a similar statistic as our satellite quenched frac- 
tion excess to show that this excess is independent of satellite 
mass. We have shown that this results extends to more ro- 
bust measures of SFR, and to all halo masses and at all radii. 
However, our physical interpretation of this statistic differs 
from theirs. In particular, this statistic only measures the 
excess fraction of satellites that quenched after infall that 
would not have quenched had they remained central galax- 
ies, because the 2 ~ central galaxy quenched fraction does 
not represent the initial conditions of satellites at the time 
of infall. We will explore this issue further in [Paper IH] 

The mass and radius dependences we see disagree with 



those of van den Bosch et al. ( 2008| b, who found only weak 
dependence of galaxy properties on halo mass and radius at 
fixed stellar mass. However, they examined primarily mean 
g — r color of satellites as a measure of star formation, which 
is problematic because dust reddening causes a severe bias- 
ing oi g — r color, particularly for low- mass galaxies, and the 



mean value is a less sensitive characterization of a highly bi- 
modal distribution. Indeed, we find similar results to theirs 
when examining mean g — r color. 

The satellite SSFR distribution dependencies we have 
detailed in this paper provide tight constrains on models for 
satellite-specific evolution. In [Paper III[ we will use a high- 
resolution cosmological simulation to obtain the detailed or- 
bital histories of satellites, which we will use to more rigor- 
ously constrain the timescales of satellite quenching, and in 
[Paper IVJ we will develop models for satellite-specific quench- 
ing via different mechanisms to test which reproduce the 
trends observed in this paper. 
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